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ABSTRACT 



A Taylor-Couette cell for the investigation of geometrical 
phase in acoustics has been constructed. The inner and outer 
cylinders are made of acrylic. The cell has an inner cylinder 
radius of 9.477 cm, radius ratio of 0.902, and cavity aspect 
ratio of 49.2. Two high performance d. c. motors can rotate 
the cylinders independently. The angular speed of the 
cylinders is monitored by a photo-interrupter. A 
piezoelectric polymer (PVDF) transducer mounted on the inner 
cylinder drives the acoustic field in the cavity and a 6.02 mm 
diameter electret microphone embedded in the wall of the inner 
cylinder acts as a receiver. Static measurements of the 
acoustic modes in the annular cavity show good agreement with 
theory, though measurements at different angular locations of 
the transducer reveal 2% nonuniformities. The effects of 
nonuniformities could be overcome by rotating the cylinders at 
high speeds. The performance under these preliminary tests 
indicate that an apparatus of this design is suitable for 
investigations of acoustics in rotating flows. 
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I. 



INTRODUCTION 



The motion of a fluid flow between two concentric, long 
cylinders undergoing differential rotation has been one of the 
most studied problems in hydrodynamics (Donnelly, 1991) , with 
however, numerous unanswered questions. Such fluid motion is 
called Taylor-Couette flow. 

The steady state axis-symmetric flow is an exact solution 
of the Navier-Stokes equation in this configuration. This 
flow is simply described as the superposition of solid body 
rotation plus the flow of a vortex line along the common axis 
of the two cylinders. When the inner cylinder is turned at a 
higher rate, such that the inner layers of the fluid are 
moving more rapidly than the outer ones, the centrifugal force 
pulls these inner layers outward which then break into cells 
and circulate. The fluid breaks into horizontal bands. Thus, 
the axis-symmetric flow is, in general, not stable. As the 
angular velocity of the inner cylinder is increased, the 
number of bands increases at first, then they become wavey. 
The velocity of rotation of this wavey pattern approaches a 
numerical constant (which is determined by the geometrical 
parameters of the system) times the angular velocity of the 
cylinder. Finally, for higher values of the angular velocity, 
or in the configuration in which the outer cylinder is 
rotating in the opposite direction, the flow becomes turbulent. 
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The complexity of behavior in this simple system is in 
principle contained in the equations of a viscous fluid. 
However, there does not exist, up to the present, a 
satisfactory method to analyze and explain the complexities 
hidden in the Navier-Stokes equation: instabilities, the onset 
of turbulent flows, and intermittency . 

Consider however, the flow not as a dynamical system, but 
as a collection of parameters of an acoustic field. In the 
geometrical acoustic approximation, the acoustic field is 
described by an equation analogous to Schrodinger' s equation 
of a charged particle, with the flow velocity playing the role 
of the magnetic vector potential. Then, it is possible to 
determine global properties of the flow by measurements of the 
acoustic phase. Modal resonance splitting and shifting are 
manifestations of global changes. 

This thesis describes the design, construction, and 
preliminary testing of a Taylor-Couette cell. The flow 
apparatus is designed with maximum flexibility in mind. 
Preliminary testing and experimentation using the cell 
determined its suitability for future research. 

The organization of this thesis is as follows. Chapter II 
presents the theory for the mode structure of the Taylor- 
Couette cell in its static and differential rotation 
configurations. Chapter III describes the design and 
construction of the cell. Chapter IV describes the 
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preliminary 
summary and 



tests and their results. Chapter V presents a 
conclusions of this work. 
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II. 



THEORY 



A. EXPECTED RESONANCE FREQUENCIES 

In order to observe the modal frequency split and shift in 
a rotating circular Couette system, we require knowledge of 
the modal structure. The following theory derives those 
frequencies . 

Consider the annular cavity formed by two concentric 
cylinders of inner and outer radii, R l and Ro, respectively, 
with top and bottom boundaries located at z= 0 and z= L z . If 
all surfaces of the cavity are perfectly rigid, so that the 
particle velocity normal to the surface is zero, then for 

linear acoustics fi • Vp=o at all boundaries, where Vp is the 



gradient of the acoustic pressure. Thus 




( 1 ) 



For a system with circular symmetry, it is best to recast 
the linear inviscid wave equation in polar coordinates 



yp + _i dp + i y P + yp _ 1 yp 

dr 2 r dr z 2 80 2 9z 2 c 2 9 1 2 ' 
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where c is the speed of sound of the medium. The technique of 
separation of variables is use to solve Eq. (2) . The acoustic 
pressure is assumed to be a product of four terms, each a 
function of only one variable, i.e. 

p(r,0,z,fc) = i?(r) 8 (0) Z(z) e Jut . (3) 



Substitution of Eq . (3) into Eq. (2) yields 



QZ— — + ® ZdR + + 

dr 2 rdr r 2 d0 2 



+ k 2 RQZ 

dz 2 



0 , 



(4) 



where k=0)/c. Multiplying each term in this equation by 1/R0Z 
and grouping terms containing z on one side results in 



d 2 i? + IdR + ld 2 Q + k2 = -ld 2 Z (5) 

Rdr 2 rRdr r 2 0d0 2 Zdz 2 

The left-hand side of this equation, a function of r and 0, 
cannot be equal to the right-hand side, a function of z alone, 
unless both functions equal some constant. If we let this 
constant be k z 2 , the right-hand side becomes 

(^ + JcJ ) z = 0 - (6) 

Repeated performance of the separation of variables technique 
applied to Eq. (4) results in the final set of equations: 
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d 2 

dz 2 



Z = 0, 



(7a) 




d6 2 



+ m 2 9 = 0 , 



(7b) 



dr 2 rdr 





R = 0, 



(7C) 



where m 2 is an arbitrary separation constant and k m 2 =k 2 -k z 2 . 
The complete general solution to Eg. (7a) is 



where A x and A 2 are arbitrary constants. The solution to Eq. 
(7b) is 



where B is an arbitrary constant and <J> is an arbitrary phase 
angle dependent on the location of the source. Eq. (7c) is 
known as Bessel's equation. Solutions to this equation are 
the transcendental functions called Bessel functions of the 
first kind J ra (k m r) and second kind Y m (k m r) of order m. 



where Ci and C 2 are arbitrary constants. Bessel functions are 
oscillating functions whose amplitudes decrease as their 



Z(z) = j^cos (k z z) + A 2 sin{k z z) , 



( 8 ) 



0(0) = Bcos(rr& + <(>) 



(9) 



R(r) = + C^Jk^) , 



( 10 ) 
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arguments increase. The second kind functions Y m (k m r), also 
known as Neumann functions, become unbound in the limit as k m r 
approaches zero. 

Applying the appropriate boundary conditions to Eq. (8) 
results in 

Z(z) = A x cos ( k z z ) , 

where 

k z = JSL {=o,i,2,... 

L z 

If the acoustic pressure p is to be a single-valued function 
of position, then p(r,0,z,t) must equal p (r, 0+2rc, z, t ) , which 
restricts the separation constant m to integral values; 
m=0, ±1, ±2, . . . 

When k m =0 and m=0, Eq. (10) implies the solution 

R(r) = C x = constant . (12) 

When k m is non-zero, the boundary conditions applied to Eq. 
(10) results in the equality 



Ym(k m Rj Ymlk a R 0 ) ' 



where the prime denotes the first derivative with respect to 
kr evaluated at the specified argument. If we let the ratio 
R 0 /Ri = 9 r then Eq. (13) can be written as 
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Ym(k m Ri) 



Y'(gk m R x ) 



( 14 ) 



The roots of Eq. (14), values of can be solved for using 
a computer for various values of m (see Appendix) . If the 
values of the m th order Bessel functions which cause Eq. (14) 
to equal zero are designated by j ran , then the allowed values 
of k m assume discrete values given by k mn = j^/R* . 

Appropriate solutions to Eq. (3) for the case k mn not equal 
to zero yield the normal modes 



Pimn (Y » 0 , Z , t) - 



imn 



Y m ( k mn r) - 

in ' mn * 






J^ k mn r ) 



( 15 ) 



x cos (/n0 + 4>) cos (k z) e JaimC , 



where 




The allowed frequencies are 



k ^ 

A-mn T •*'z 

(= 0 , 1 , 2 ... 

m = 0,1,2... 

n = 0,1,2... 
m = n * 0 . 

( 16 ) 

ius quantized by the relationship 
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( 17 ) 



~\mn 




Each eigenfunction given by Eq. (15) has its own 
characteristic eigenf requency, denoted by Eq. (17) , and can 
be specified by the ordered integers (t,m,n) . The value of t 
determines the number of nodal surfaces in pressure 
perpendicular to the z axis. The value of m specifies the 
number of radial pressure nodes. The value of n gives the 
number of nodal surfaces of constant r (cylinders) surrounding 
the z axis and concentric with it. 

Solutions to Eq. (3) when k mn equals zero yield the normal 
modes 

P t ,o,o (r,6,z,t) = P, 0 O cos (k z z) e jat , 



where 



* 2 - = k 2 z 

c 2 

k z = - T - « - 0,1,2, . . . 



( 18 ) 



The allowed frequencies are 



-i,o,o 




corresponding to pure z modes. 



( 19 ) 
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B. MODE STRUCTURE FOR DIFFERENTIAL ROTATION 



Axis-symmetric flows in a Taylor-Couette cell are derived 
for the annular cavity between two infinitely long cylinders. 
However, the expected resonance frequencies derived assumed 
finite length cylinders. A justification for assuming such 
velocity field has been given previously by Andereck et al., 
1983. These researchers noted that if the aspect ratio T, the 
ratio of the fluid height to the gap between the cylinders, at 
the region of measurement is properly chosen, then the 
boundary conditions existing at the top and bottom of the 
cylinders exerted no observable difference in the flow state 
properties. In other words, boundary conditions at the top 
and bottom of the cylinders can be neglected if measurements 
are made far removed from these boundaries. Many of the 
measurements cited in the above study were conducted with an 
aspect ratio of 30. However, other measurements were made 
with 20<r<47. Our initial testing is conducted with T=24.6. 
Therefore, we can assume that our region of interest has an 
azimuthally symmetric velocity field. 

Consider waves with high azimuthal numbers m. For these 
waves the velocity field induced by the cylinders' 
differential rotation varies over distances large compared 
with the wavelength. Within this approximation, the local 
dispersion relation in the moving medium with velocity flow 

u(r) is given by 
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0) = a) ( k) + k ‘u ( r) , 



( 20 ) 



where u ( k ) is the dispersion relation of a wave with wave 



vector £ in the stationary medium. A particular realization 



of the flow velocity is given by the stable azimuthal 
configuration with a stationary inner cylinder. For flow 
velocities I u I c c, Eq. (20) can be written to lowest order as 

0) = a>(|ic+ (21) 



The correspondence between 0) and the Hamiltonian establishes 
the analogy between the magnetic vector potential A and the 
flow velocity, provided 



Q 

h 




( 22 ) 



(Berry et al . , 1980; Larraza, 1989) Making the natural 

identification between the wave vector k and the operator 



-iV, Eq. (21) is analogous to Schrodinger ' s equation for a 



charged particle in a vector potential A. The geometry of 

Taylor-Couette flow for the azimuthally symmetric case yields 
the eigenf requency spectrum 
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m + 



( 23 ) 



ca 

c 




MQg 0 

c 



(Larraza, 1989) , valid for rotation of the outer cylinder 
only. If the rotational speed £2=0, the states m*0 are doubly 
degenerate. For a non-zero value of the rotational speed the 
degeneracy is broken. Physically, this means that the waves 
travelling in the same direction as the average flow have an 
additional contribution to the sound speed, whereas waves 
travelling in opposite direction to the flow have an effective 
sound speed which is less than the static value. According to 
Eq. (23) , there is an asymmetric splitting of the degenerate 
states about the zero-angular speed eigenfrequency . To see 
this, notice that Eq. (23) is the first order expansion of 



c 




( 24 ) 



„ l/n|£2.R„ n . , , „ 

For 2 — ! £ = _ with n an integer, degeneracy occurs. For n=l, 

C 2 

the states m=l and m=-2 are degenerate with frequency 
3C 



CO =- 



2R r 



, while the state m=-l is not degenerate. Physically, 



these states are reached by increasing the rotational speed £2 
and correspond to high azimuthal numbers m. 
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Ill . 



TAYLOR-COUETTE FLOW APPARATUS 



A. MECHANICAL DESIGN 

With flexibility in mind, the cell is designed to allow 
the inner and outer cylinders to be rotated independently. 
Also, rigid boundaries at the top and bottom of the cylinders 
allow for three types of boundary motion: stationary, rotation 
at the angular speed of the inner cylinder, and rotation at 
the angular speed of the outer cylinder. The system is 
designed to be vertically mounted to eliminate flexural strain 
resulting from the weight of the components if mounted in a 
horizontal configuration. 

1. Cylinders 

The inner and outer cylinders are made of acrylic. 
The acrylic is fabricated using extruded, vice cast, methods, 
which introduced nonuniformity effects which will be described 
later. The inner cylinder radius^, is 9.477 cm and the 
outer cylinder radius, Ro, is 10.510 cm. The resulting radius 
ratio is 0.902, which is near the value used in many previous 
studies of Taylor-Couette flow. (DiPrima and Swinney, 1981) 
The solid boundaries used to enclose the annulus create a 
resonant cavity height, L*, of 50.838 cm, resulting in a 
cavity aspect ratio of 49.2. 
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2 . Support Structure and Bearings 



With the exception of the cylinders, drive motors and 
their associated drive belts, the entire assembly, including 
support structures, are made of aluminum. The inner cylinder 
is supported by a length of solid round stock, with a diameter 
of 5.08 cm, which also provides axial alignment for the 
bearings . 

Radial bearings are used to support rotating 
components since radial rigidity is of paramount importance. 
The type of radial bearings chosen are able to handle compound 
loads. Thus, while providing rigidity, they are also able to 
support load resulting from the weight of the assembly. 

3 . Upper and Lower Boundaries 

Aluminum rings provide an acoustically rigid boundary 
condition. The rings can either be stationary, or rotated 
with the same angular speed as the cylinder to which they are 
attached. The rings are machined to fit flush against the 
cylinder they are connected to while providing as narrow a gap 
as possible near the other cylinder. While a tight cavity is 
desirable to set up standing waves, some gap between the rings 
and the cylinder walls is necessary to eliminate frictional 
rubbing which results in torque disturbance and noise pick-up 
on the receiver. The average gap between ring and cylinder 
wall is 0.01 in. 
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4 . Motors and Drive Belts 



Two high performance, fractional horsepower, disc- 
armature, dc motors are used to rotate the cylinders. These 
motors, PMI U-series model #U12M4H, are designed for excellent 
speed control and constant torque delivery. The motors are 
connected to the geared cylinders via a stainless steel cable 
reinforced, polyurethane, toothed drive belt. 

5. Miscellaneous 

The entire assembly is mounted on a 2 cm thick 
aluminum baseplate measuring 50 cm by 62 cm. Two vertical 
aluminum cylinders support a horizontal crossbar which 
provides axial alignment of the components at the top of the 
assembly. Figures 1, 2, and 3 show the construction of the 
apparatus as well as close-ups of the drive assembly and top 
cylinder mounting. 
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Figure 1 Taylor-Couette cell flow apparatus 
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Figure 2 Expanded view of the drive assembly 
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Figure 3 Expanded view of upper mounting 
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B. ELECTRICAL DESIGN 



1 . Transduction 

Several constraints, mechanical as well as 
theoretical, led us to choose a piezoelectric polymer film as 
the source of cavity excitation. The polymer is 
polyvinylidene fluoride (PVDF) . PVDF is sold commercially in 
sheet form, resembling a piece of thick plastic wrap. For our 
application, a 4 0|lm thick strip was cut 2.0 cm wide and 50.8 
cm long to match the length of the cavity. For a constant 
driving voltage, PVDF acts as a constant displacement source. 

Our ultimate goal is to investigate changes in 
azimuthal mode structure due to the rotating flow. 
Consequently, it is not necessary to excite radial and 
vertical modes in order to observe the frequency split or 
shift. Choice of excitation frequency and source placement 
can limit excitation of these modes. 

Since the PVDF strip is a constant displacement 
source, its width has to be less than one-half wavelength for 
the range of frequencies of interest. The upper frequency 
limit for preliminary testing was chosen to be 8kHz. Similar 
reasoning as that used for strip width was also applied to 
strip placement. With the strip extending across the length 
of the cavity, those vertical modes requiring a node at the 
mid-height of the cavity would be weakly excited, if at all. 
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2 . Receiver 



A 6.02 mm diameter electret microphone is embedded in 
the wall of the inner cylinder so that its face is flush with 
the surface of the cylinder. A silicone gel is placed around 
the microphone to provide a tight seal. The desire to 
minimize the number of penetrations through the cylinders 
during preliminary testing led us to use only one microphone 
placed at the vertical center of the cavity, 25.4 cm from 
either vertical boundary. However, the use of only one 
microphone and its placement necessarily limits the 
interpretation of the results of the preliminary testing as 
will be discussed in Chapter IV. 

3 . Angular Speed 

The design uses a photo-interrupter to measure the 
angular speed of the outer cylinder. A small tab, mounted at 
the top of the outer cylinder, gates the output of the photo- 
interrupter once each revolution. The interrupter signal 
output is fed to an HP 5384A frequency counter which counts 
the number of events occurring within a .1, 1, or 10 second 
window. Digital readout could either be in hertz or the 
average period of revolution. 
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IV. 



RESULTS OF PRELIMINARY TESTS 



To assess the capabilities and limitations of our initial 
design, we performed two preliminary tests. The purpose of 
these tests was to verify our theoretical analysis of the cell 
and to determine the magnitude of errors caused by 
nonuniformities . 

With the inner and outer cylinders fixed, the output of an 
HP 3562A dynamic signal analyzer was used to excite the 
cavity. A 5 Vrms periodic chirp source output, a fast sine 
sweep over the selected frequency span that repeats with the 
same period as the time record, was used. The output of the 
embedded electret microphone was fed back to the signal 
analyzer and its power spectrum obtained. 

Frequencies corresponding to clearly discernable peaks in 
the power spectrum were recorded. The power spectrum was 
analyzed from d.c. to approximately 3.6 kHz in 400 Hz windows, 
having 0.5 Hz resolution. Within this frequency span, 
approximately sixty theoretical eigenf requencies are possible. 
To aid in the resonant frequency cataloguing, an HP 3314A 
function generator was used in conjunction with the dynamic 
signal generator to manually sweep through a resonance and 
observe the phase change on an oscilloscope. 

The results of these measurements are shown in Figure 4, 
which is a graph of observed resonance frequencies versus the 



21 



corresponding theoretical frequencies. The solid line has a 
slope of one. Ideally, the data points should fall on this 
line. The degree to which they do not is shown in Figure 6, 
which shows the magnitude of percent error. Generally, the 
results are quite good. The agreement is typically within 1%. 

A number of modes were easily observable indicating the 
gap between the boundary rings and the cylinder walls did not 
kill the mode structure. The mode structures of interest, 
those with purely azimuthal dependence, were often masked by 
much stronger, closely spaced vertical modes. Vertical mode 
structures were not eliminated to the degree we had 
anticipated. This maybe due in part to the way the PVDF is 
attached. The strip is actually composed of two shorter 
strips, joined at the midpoint of the cavity. The joint 
construction may result in a smaller displacement of the PVDF 
at the joint as compared to the rest of the strip, therefore, 
allowing vertical modes to be excited. Further work in this 
area is needed. 

The use of one microphone, located at the center of the 
cavity, is not sufficient to distiguish between closely spaced 
modes. A moveable probe microphone might be better suited for 
determining mode structure. 

Lastly, we did not observe any purely azimuthal modes and 
only about half of the modes represented in Figure 4 have any 
azimuthal component. This again points to the need for 
continued work on the method of transduction. 
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Figure 4 Observable resonance peaks from dc to 
3.6kHz 




Figure 5 Agreement between observed resonance and 



theoretical resonance 
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While making the measurements discussed above, we noticed 
that the resonance frequency of a mode depended on the 
relative angular orientation of the inner and outer cylinders. 
We suspect that this is due to nonuniformities in the cylinder 
walls, ring boundaries, and alignment. The magnitude of this 
effect is illustrated in Figure 6. The data show the result 
of measurements of the resonance frequency at six equally 
spaced relative orientations. In Figure 6, the percentage 
variation of the resonance frequencies is plotted as a 
function of the mean resonance frequency. In all cases, the 
percentage variation is less than approximately 1.5%. This 
could account for much of the error shown in Figure 5. 
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Percent Variation 




Figure 6 Variation of resonance through 360 degrees of 
rotation 
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V. 



SUMMARY AND CONCLUSIONS 



This thesis reports the design and construction of a 
Taylor-Couette cell. The apparatus uses acoustic transduction 
as a means to characterize global properties of the flow. 
Preliminary tests show that the static mode structure differs 
slightly from theoretical predictions, mainly because 
nonuniformities in the cavity. The nonuniformities were 
qualitatively observed by measuring the resonant frequencies 
for different angular positions of the transducer. These 
preliminary tests provide qualitative limitations imposed by 
the apparatus for the applicability of the theory in Chapter 
II. Because of nonuniformities, certain degeneracies are 
asymmetrically removed, mocking the effect of an arbitrary 
flow. Thus, quantitative knowledge of the effect due to 
nonuniformities is essential in order to determine the minimum 
differential rotation for which the effects due to the flow 
become important. 

The acoustic "visualization" reported in this thesis 
constitutes a novel concept, but further development of the 
technique is required. As opposed to standard visualization 
methods, in principle acoustic transduction does not limit the 
Reynold's numbers for which the technique is useful. But most 
important, because the acoustic phase can provide global 
properties of an arbitrary flow, the investigations proposed 
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in this thesis may lay the foundations for a first-principles 



phenomenological theory of turbulence that relies 
complete set of global variables for its description. 



on a 
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28.8 


3,1 


28.0 


0,2 


57.6 


1,2 


57 . 6 
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MODES OF A CYLINDRICAL ANNULUS 



m , n \ 1 


0 


• 


2 


4 


4 


5 


6 


0,0 


0 


339 


678 


1016 


1355 


1694 


2033 


1,0 


549 


645 


872 


1155 


1462 


1781 


2106 


1 2,0 


1098 


1149 


1290 


1496 


1744 


2019 


2310 


4,0 


1646 


1681 


1780 


1935 


2132 


2362 


2616 


4,0 


2195 


2221 


2298 


2419 


2580 


2773 


2992 


5,0 


2744 


2765 


2826 


2926 


3060 


3224 


3415 


4,0 


3293 


3310 


3362 


3446 


3561 


3703 


3870 


7,0 


3842 


3857 


3901 


3974 


4074 


4199 


4346 


4,0 


4390 


4403 


4442 


4506 


4593 


4706 


4839 


3,1 


16661 


16664 


16674 


16692 


16716 


16746 


16784 


1,1 


16661 


16664 


16674 


16692 


16716 


16746 


16784 


2,1 


16716 


16722 


16732 


16749 


16773 


16804 


16842 


3,1 


1 6776 


16780 












0,2 


33321 


33323 


33328 










1,2 


33321 


33323 


33328 











Cavity Height 



Inner Radius 



Outer Radius 



Ro/Rt 



50.838 cm 



9.477 cm 



10.510 cm 



1 .109 
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